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This  study  was  undertaken  as  part  of  an  investigation  of  navi¬ 
gation  within  the  solar  system  by  optical  means.  The  objective  of 
the  investigation  is  to  evaluate  the  suitability  of  various  physical 
phenomena  as  sources  of  navigational  information  and  to  estimate  the 
accuracy  of  navigational  information  obtained  by  various  techniques. 
The  work  was  supported  under  Contract  No.  AF  33(616) -7413  by  the 
Navigational  and  Guidance  Laboratory,  Wright  Air  Development  Division, 
Air  Research  and  Development  Connand,  United  States  Air  Force. 

The  present  report  is  one  of  a  series  of  three  entitled: 

SATELLITE  NAVIGATION  BY  TERRESTRIAL  OCCULTATIONS  OF  STARS 

I.  General  Considerations  Neglecting  Atmospheric  Refraction  and 
Extinction  by  Ali  M.  Naqvi  (Geophysics  Corporation  of 
America  Technical  Report  62-18-A,  October  1962). 

II.  Considerations  Relating  to  Refraction  and  Extinction  by 
Ali  M.  Naqvi  (Geophysics  Corporation  of  America  Technical 
Report.  62-21-A,  October  1962). 

III.  Interference  Due  to  Brightness  of  Earth's  Atmosphere  by 

Rollin  C.  Jones  and  Ali  M.  Naqvi  (Geophysics  Corporation  of 
America  Technical  Report  62-22-A,  October  1962). 
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ABSTRACT 


The  condition  for  terrestrial  occultation  of  a  star  as  seen  from 
a  satellite  is  derived  in  terns  of  the  orbital  elements  of  the  sat¬ 
ellite  and  the  equatorial  coordinates  cf  the  star.  An  occultation 
equation  relating  the  time  of  occultation  to  the  orbital  elements  is 
also  derived.  Relevant  stellar  properties  are  discussed.  The 
frequency  of  occultations  as  a  function  of  perigee  distance  and 
eccentricity  of  the  orbit  is  obtained.  A  general  discussion  of  the 
minimum  number  of  stars  required  for  satellite  navigation,  and  the 
corresponding  limiting  magnitude  of  the  satellite-borne  telescope  is 
given.  Except  for  references  to  other  reports  of  this  series  where 
atmospheric  effects  are  considered,  all  results  are  derived  on  the 
assumption  that  the  earth  is  a  perfect  sphere  and  has  no  atmosphere. 
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SECTION  1 

INTRODUCTION 

One  g£  the  classic*!  method*  for  measuring  the  longitude  of  a 
place  on  the  earth's  surface  is  based  upon  the  observations  of  star 
occultations  by  the  moon  (See  Chauvenet,  1891,  Chap.  10).  The  method 
assumes  a  knowledge  of  the  lunar  orbit.  Alternatively  if  the  longi¬ 
tude  of  the  observer  could  be  determined  by  some  other  method,  the 
star  occultations  could  be  used  to  determine  the  lunar  orbit,  or 
rather  to  apply  small  corrections  to  an  assumed  orbit.  The  similarity 
between  this  procedure  and  a  proposed  method  for  determining  the  orbit 
of  an  earth-satellite,  based  upon  the  measurements  of  the  time  of 
occultations  of  stars  by  the  earth,  is  obvious. 

A  parallel  beam  of  light  from  a  star  falls  on  the  earth,  which 
caste  a  cylinderical  shadow  whose  axis  is  along  the  line  Joining  the 
star  and  the  center  of  the  earth.  As  seen  from  the  satellite  an 
occultation  of  this  star  by  the  earth  will  occur,  if  and  only  if, 
some  portion  of  the  satellite's  orbit  passes  through  the  above  mention¬ 
ed  cylinderical  shadow. 

The  earth  is,  in  fact,  not  a  perfect  sphere  and  therefore  its 
shadow  is  not  a  perfect  cylinder.  However,  correction  for  Known 
departures  from  sphericity  can  be  applied.  To  the  extent  that  there 
are  bound  to  be  some  errors  in  the  measurements  of  the  shape  of  the 
earth,  the  exact  shape  of  the  shadow  would  also  be  uncertain.  Then 
there  are  mountains  on  the  earth's  surface  some  of  which  extend  to 
almost.  10  km  above  the  earth's  surface.  Unless  their  exact  locations 
and  height  are  known,  some  small  errors  can  arise  due  to  these  irreg¬ 
ularities. 
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However,  the  presence  of  clouds  in  the  earth's  atmosphere,  which 
can  completely  absorb  the  star  light,  and  of  which  the  height  may 
vary  in  a  totally  unpredictable  manner,  makes  the  observations  of 
occupations  of  stars  by  the  earth's  solid  surface  unreliable  for 
position  and  orbit  measurements.  In  the  second  report  of  this  series 
(Naqvi,  1962;  hereafter  referred  to  as  Report  II),  it  is  shown  that 
criteria  other  than  a  geometrical  occupation  behind  the  opaque  earth 
can  be  used  for  satellite  navigation.  In  this  report  we  shall  assume 
that  the  earth  is  a  perfect  sphere  with  no  atmosphere  (and  no  clouds) . 
Section  2  defines  the  various  orbital  elements.  In  Section  3  a  condi f 
tion  for  the  occupation  of  a  given  star  by  a  satellite  of  a  given 
orbit  is  derived.  Usually  an  approximate  orbit  of  the  satellite  is 
predetermined.  This  condition  is  then  useful  since  it  allows  us  to 
know  in  advance  of  the  launching  of  uhe  satellite  exactly  which  stars, 
(and  how  many),  would  be  occulted.  In  Section  4,  the  occupation 
equation  which  relates  the  time  of  occupation  to  the  orbital  elements 
is  derived.  A  resume  of  relevant  stellar  properties  is  given  in  Sec¬ 
tion  5,  and  finally,  the  expected  frequency  of  star  occultations  and 
the  limiting  magnitudes  of  the  stars  to  be  observed  are  discussed  in 
Section  6. 

In  Report  II  of  this  series,  the  effects  of  refraction  and 
extinction  are  Investigated  and  in  Report  III  a  brief  discussion  of 
the  interference  due  to  the  brightness  c£  the  .earth's  atmosphere 
is  given. 
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SECTION  2 


THE  ORBITAL  ELEMENTS 

In  this  section  we  define  the  six  fundaments!  orbital  elements 
and  a  few  others  related  parameters.  For  more  complete  details  the 
reader  may  consult  any  standard  text  on  celestial  mechanics  or  orbit 
theory  (for  example  Moulton,  1914).  For  an  elliptical  orbit  six 
orbital  elements  and  time  completely  define  the  position  of  a  body  in 
space.  (Figure  la).  These  are: 

a  •  semi-major  axis,  which  defines  the  size  of  the  orbit. 

e  -  eccentricity,  which  defines  the  shape  or  oblateness  of  the 
orbit. 

T  -  time  of  perihelion  passage  which  defines  the  position  of  the 
body  in  its  orbit  at  a  given  time, 
i  -  inclination  of  the  orbital  plane  to  the  fundamental  plane 

(xy  plane  in  Figure  la).  Note  that  in  the  usual  orbit  theory 
the  ecliptic  is  taken  to  be  the  fundamental  plane,  but  in 
considering  the  satellite  orbits,  the  earth's  equatorial  plane 
is  found  more  convenient. 

fi  -  longitude  (if  the  ecliptic  coordinate  system  Is  used)  or  right 
ascension  (if  the  equatorial  coordinate  system  is  used)  of  the 
ascending  node.  The  ascending  node  is  defined  as  the  point 
of  intersection  of  the  orbit  and  the  fundamental  plane,  when 
the  body  is  moving  from  negative  to  positive  Z  direction  and 
i  and  (1  together  define  the  orientation  of  the  orbital  plane 
with  reference  to  a  given  coordinate  system. 
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Figure  la.  The  orbital  elements  in  the  equatorial  coordinate  system, 
(xy-plane  is  the  plane  of  the  earth's  equator). 


Figure  lb.  The  orbital  elements  in  the  orbital  plane. 
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u  -  longitude  of  the  perigee  measured  fro*  the  ascending  node 

along  the  orbital  plane,  which  defines  the  orientation  of  the 
orbit  in  its  plane. 


There  are  several  other  orbital  parameters  which  are  functions 
of  the  above  mentioned  elements,  and  which  can  therefore  be  used  equal¬ 
ly  well  in  their  place.  He  shall  have  occasion  to  use  two  such  para¬ 
meters,  which  are: 

P  -  period  of  the  body  which  is  related  to  the  semi-major  axis  a 
in  accordance  with  Kepler's  third  law; 


3 

a 


2  ■> 

4 «  1  4* 

— —  ■-'■■■  ^ 

G(m  +  m1)  Ga 


constant 


(2.1) 


The  approximation  in  Equation  (2.1)  holds  provided  the  kiss 
of  the  smaller  body  m'  is  negligible  compared  to  the  mass  of 
the  larger  body  m.  G  is  the  gravitational  constant.  A 
convenient  way  to  use  Kepler's  law  is  to  write  it  as  follows: 

*  3/2 

P  -  P0  (  /*Q)  (2.2) 

where  PQ  is  the  period  of  an  imaginary  satellite  whose  semi- 
major  axis  is  iq.  If  ao  ■  6378  km,  the  equatorial  radius  of 
the  earth,  we  have 

PQ  "  5069  sec  ■  84.48  min  ■  1.408  hours  (2.3) 

v  -  true  anomaly  is  the  polar  angle  (Figure  lb)  measured  from  the 
perigee  in  the  forward  direction  of  motion  along  the  orbit. 

It  is  evident  that  v  is  independent  of  Che  elements  which 
define  the  orientation  of  the  orbital  plane  in  space,  or  the 
orbital  orientation  in  this  plane,  namely  i,  0,  and  u,  and  is 
a  function  of  a  (or  P),  e,  T  and  the  time  t.  For  derivation 
of  the  following  results  a  reference  may  be  made  to  Moulton 
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(Chap.  5: 


1914). 


Cap 


v  ll  +  e; 

2~  F"*i 


4  % 


tan# 

2 


(2.4) 


E  -e  sin  E  *  M  (known  as  Kepler's  equation)  (2-5) 


"  ■ &  M  «.»> 

When  the  eccentricity  e  is  sms  11  (as  is  the  case  for  planetary 
orbits,  and  for  some  satellite  orbits)  a  series  expansion  in  e  can 
be  obtained 


v  -  M  +  2e  sin  H  +  |  e2  sin  2M 

3 

+  ^  (13  sin  3M  -  3  sin  K) 

4 

+  |jr  (103  sin  4H  -  44  sin  2H) 

5 

+  (1G97  sin  5M  -  645  sin  3H  +  50  sin  M) 

6 

+  §50  (1223  sin  6M  -  302  sin  4M  +  85  sin  2M) 

+ .  (2.7) 

The  polar  equation  of  an  ellipse  in  its  simplest  for"  is 
=  a(l  -  e2> 

r  1  +  e  cos  v  (2.8) 
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SECTIOK  3 

COREmORS  FOR  OOCCLTATIOR 

It  can  be  seen  from  elementary  geometrical  consideration  that  if 
the  earth-star  vector  is  in  the  plane  of  the  satellite  orbit,  an 
eccultstics  of  this  star  ss:  occur  (See  Figure  2}  during  an  orbital 
revolution.  On  the  contrary,  if  the  earth-star  vector  is  perpendicu¬ 
lar  to  the  orbital  plane  an  occultation  can  not  occur.  The  general 
condition  is  derived  in  this  section. 


There  is  another  way  to  look  at  the  occultation  geometry  froa  the 
point  of  view  of  an  observer  on  the  satellite,  which  is  depicted  in 
Figure  3.  A  cone  is  drawn  just  enveloping  the  earth  with  the  satellite 
at  its  apex  and  the  line  joining  the  satellite  to  the  earth’s  center 
as  its  axis.  The  semi-angle  of  the ’Shadow  cone'*  is 


0  *  sin 


-1 


a 

o 

r 


(3.1) 


where  a  is  the  earth's  radius  and  r,  the  distance  from  the  earth’s 
o 

center  to  the  satellite,  r  is  sometimes  expressed  in  terms  of  h, 
the  height  above  the  earth’s  surface,  as 


r  *  aQ  +  h  (3.2) 

A  star  whose  direction,  as  viewed  from  the  satellite,  lies  on  the 
forward  half  of  the  instantaneous  surface  of  the  cone  is  said  to  be 
at  the  point  of  ingress,  whereas  a  star  whose  direction  lies  on  the 
other  half  of  the  instantaneous  surface  of  the  cone  is  said  to  be  at 
the  point  of  egress. 
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Figure  3.  Occultation 


geometry 


As  Che  satellite  revolves  -round  the  earth,  the  "shadow  cone", 
which  can  be  imagined  to  be  attached  to  the  satellite,  also  revolves. 
For  an  elliptical  orLit  the  radius  vector  r  is  given  by  Equation  (2.8). 
Variations  in  r  cause  changes  in  the  seai-angle  of  the  shadow  cone  p, 
which  as  a  function  ot  the  true  anomaly  v,  becomes. 


.  -1  o  1  +  e  cos  v 
e«sin  - - - - — 

1  -  e 

For  any  orbit  the  maximum  and  minimum  values  of  p  are 


P(maj)  *  sin 


-1 


a(l  -  e) 

p(aih)  «  sin'1 


(3.3) 


(3-4) 

(3.5) 


These  are  functions  of  only  two  of  the  orbital  coordinates,  namely  a 
and  e.  It  is  evident  that  the  condition  for  occultation  is  that  the 
line  joining  the  star  and  the  satellite  must  pass  through  the  shadow 
cone  at  some  point  along  the  satellite  orbit. 


Consider  a  celestial  sphere  (Figure  4)  with  the  earth  at  its 
center  and  the  plane  of  the  satellite's  orbit,  BSPQS',  as  its  funda¬ 
mental  plane.  Let  F  denote  the  perigee  and  S  the  projection  on  the 
celestial  sphere  of  the  instantaneous  position  of  the  satellite  in  its 
orbit.  The  "shadow  cone",  described  above  would  intersect  the  celes¬ 
tial  sphere  in  the  shaded  circle  whose  center  S'  is  *  radians  away 
from  S.  The  angular  radius  of  this  circle  is  defined  by  Equation 
(3.3). 


Consider  a  star  X  whose  latitude  is  ijr  and  whose  longitude  $ 
(measured  from  perigee  P  along  the  orbital  plane  in  the  direction  of 
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Che  satellite's  Motion)  is  the  same  as  the  instantaneous  longitude  of 
S'.  This  system  of  coordinates  will  be  called  the  "orbital"  coordi¬ 
nate  system.  The  condition  that  an  occultation  of  this  star,  as  seen 
by  a  satellite  S  (of  given  orbit),  would  occur  is  that 


i  sin 


_1  *e(i  -  e  cos  ♦) 
e(l  -  e2) 


(3.6) 


Since  stellar  positions  are  usually  given  in  right-  ascension 
and  declination,  some  relations  between  these  coordinates  and  the 
orbital  latitude  and  longitude  are  required.  To  derive  these  rela¬ 
tions,  consider  the  celestial  sphere  drawn  in  Figure  5,  in  which  the 
plane  of  the  earth's  equator  is  the  fundamental  plane.  This  differs 
from  the  usual  drawings  of  orbits  of  bodies  in  the  solar  system  which 
have  the  ecliptic  as  the  fundamental  plane  of  the  celestial  sphere. 

The  following  results  can  be  obtained  very  easily  by  standard 
methods  of  spherical  trigonometry. 

sin  \|r  *  -sin  i  cos  6  sin  (a  -  ft)  +  cos  i  sin  5  (3.7) 

cos  <Jr  sin  (u  +  0)  *  cos  i  cos  5  sin  (a  -  (1) 

+  sin  i  sin  &  (3.3) 

cos  Jr  cos  (u  +  <>)  -  cos  (a  -  0)  cos  6  (3.9) 

Equations  (3.8)  and  (3.9)  can  be  further  simplified  to  yield 

cos  $  sin  ❖  *  cos  i  cos  6  sin  (a  -  fi)  cos  u  +  sin  i  sin  5  cos  u 

-  cos  (a  -  n)  cos  6  sin  u  (3.10) 
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Figure  5.  The  equatorial  and  the  "orbital"  coordinate 

sys tows . 
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COS  i  cos  ♦ 


■  cos  i  cos  5  sin  (a  -  n) 
+  cos  (CC  -  n)  cos  8  cos 


sin  u  +  sin  i  sin  8  sin  u 
u  (3.11) 


The  condition  for  occultstion  becomes 


—  (1  -  e^)  (cos  i  sin  8  -  sin  (ce  -  fl)  sin  i  cos  8) 
o 


■¥  e  cos 


-1  sin  1  sin  8  +  sin  (Ot  -  fl)  cos  i  cos  8 
,n  cos  (a  -  $1)  cos  6 


<  1 


(3.12) 
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SECTION  4 


DERIVATION  OF  ORBITAL  ELEMENTS 

Consider  the  celestial  sphere  drawn  in  Figure  6.  The  instantan¬ 
eous  position  of  the  satellite  is  marked  at  S  and  the  projection  of 
the  earth's  "shadow"  on  the  celestial  sphere  is  again  depicted  by 
the  shaded  circle  whose  center  S'  is  n  radians  from  S.  The  true 
anomaly  of  the  satellite  v  equals  PS.  The  star  X,  whose  orbital  lat¬ 
itude  XX'  *  jr  and  whose  orbital  longitude  PBX'  *  ®  is  at  the  point  of 
ingress. 

It  is  easy  to  see  that  XS'  («  p)  is  the  instantaneous  angular 
radius  of  the  earth's  shadow,  given  by  Equation  (3.3),  and  S'X*  * 

(c  -  «  -  v).  Now  consider  the  spherical  triangle  XX 'S'  whose  angle 
XX 'S'  is  a  right  angle.  Solving  this  triangle  we  obtain 

cos  p  ■  cos  y  cos  (4  -  n  -  v)  =  -  cos  y  cos  (v  -  «■) 

(A.l) 

We  now  use  Equations  (3.10)  and  (3.11)  to  change  from  the  orbital  to 
the  equatorial  coordinates  of  the  star,  and  get  * 

f  a  (l  +  e  cos  v)\  2  ^ 

-  cos  p  *  -II-  - » - ) 

LI  a(1  '  8  >  / 

■  cos  i  cos  6  sin  (a  -  fi)sin  (v  +  u>)  +  sin  i 

sin  6  sin  (v  +  w)  +  cos  &  cos  (a  -  fi)  cos  (v  +  •■>) 

(A.  2) 

The  observations  consist  of  the  measurement  of  the  time  of  ingress 
(or  egress)  of  several  stars.  Six  such  time  measurements  will  yield 
six  equations  from  which  the  six  orbital  elements  can  be  calculated. 

*  This  equation  has  also  oeen  derived  by  scientists  at  Control 
Data  Corporation  (Private  communication  from  Robert  Lillestrand). 


SECTION  5 


RESUME^ OF  STELLAR  PROPERTIES 

5.1  SPECTRAL  CHARACTERISTICS  AND  TEMPERATURES 

The  spectral  classification  of  stars  divides  them  into  several 
groups  labelled  B,  A,  F,  G,  K  and  M.  There  are  a  few  other  groups 
also,  labelled  W,  0,  Q_,  R,  N  and  S,  but  since  they  consist  of  only 
about  1%  of  all  the  apparently  bright  stars,  they  will  not  be  consid¬ 
ered  here.  Each  spectral  class  is  further  divided  into  10  subclasses, 
labelled  by  numbers  from  zero  to  9. 

In  addition  the  stars  are  divided  into  groups  according  to 
their  intrinsic  (as  opposed  to  apparent)  luminosity.  The  most  impor¬ 
tant  luminosity  classes  are  the  supergiants,  the  giants  and  the  main 
sequence  stars,  although  the  following  more  complete  classification 
is  frequently  used. 

I  Supergiants  (including  c  stars) 

II  Bright  giants 

III  Giants 

IV  Sub-giants 

V  Main  sequence  (also  called  dwarfs) 

Sd  Sub -dwarfs 

w  White  dwarfs 

Except  for  the  presence  of  a  few  absorption  and  emission  lines, 
the  spectral  intensity  distribution  of  a  star  closely  approximates  the 
Planck  law.  Table  1  (taken  from  Allen,  1957)  gives  the  effective 
temperatures  of  stars  of  various  spectral  and  luminosity  classes.  The 
effective  temperature  T  is  defined  in  terms  of  the  luminosity  L 

c  * 
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TABLE  i 


EFFECTIVE  TEMPERATURES  AND  COLOR  INDICES  FOR  VARIOUS  SPECTRAL  AND 
LUMINOSITY  CLASSES  (taken  from  Allen,  1957) 
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I  *  ll1  t 


:  J 

C 

in 

\+f 

(total  energy  output  per  second)  as  follows, 

2  4 

L  *  4«R  T  c 

e 

where  R  is  the  radius  of  the  star  and  a  is  the  Stefan-Boltzmenn  con¬ 
stant.  There  are  several  different  definitions  of  stellar  temp¬ 
eratures  used  by  astrophysicists  (color  temperature,  brightness 
temperature,  etc.)  but  the  effective  teaperature  is  physically  the 
most  significant,  and  it  is  a  parameter  of  great  theoretical  useful¬ 
ness  (See,  for  example,  Aller  1953,  Chap.  6). 

From  the  wide  range  of  stellar  temperatures  shown  in  Table  1,  it 
is  evident  that  the  intensity  distribution  of  different  stars  in 
S-  different  spectral  regions  would  vary  very  much.  This  is  shown  in 

Figure  7  (taken  from  a  recent  paper  by  Ramsey,  1962  *),  where  the 

j  energy  received  at  the  top  of  the  earth's  atmosphere  from  some  of  the 

V 

l  brightest  stars  (of  different  spectral  and  luminosity  classifications) 

is  plotted  against  the  wave  length.  Thus  it  can  be  seen  that  a  red 
star  which  is  visually  bright,  for  example  Betelgeux,  has  very  little 
energy  at  the  short  wave  length  region,  say  below  5000  A. 

Some  relevant  data  concerning  the  stars  for  which  curves  of  Figure 
7  are  drawn  is  listed  in  Table  2. 

The  following  stars  are  omitted  from  the  figure  of  Ramsey: 

R  Hydrae  because  of  its  low  brightness 
a  Crucis  and  Altair  because  the  temperatures  used  by 

Ramsey  are  in  error.  Their  effective  surface  tempera¬ 
tures  are  21,000*  K  and  8,600°  K  respectively,  whereas 
Ramsey  used  2,81'°  K  and  7,500°  K. 

Antare8  because  Ramsey  labelled  two  different  curves  for 
this  star.  (It  appears  that  one  of  these  curves  is  for 
Antares  and  the  other  for  Aldebtran.) 
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Figure  7.  Spectral  irradianco  of  selected  bright:  »C«r»  at  Che  Cop  of  Che 
atmosphere  (Redrawn  from  Ramsay,  1962). 


5.2  ucTEjcsrr-r  chmacieexstics 


The  range  la  the  apparent  brightnesses,  or  ^parent  magnitudes 
of  stars  (due  in  part  to  range  in  intrinsic  brightness  and  in  part 
to  that  in  distances)  Is  enormously  great.  The  naked  eye  can,  under 
optima*  conditions,  see  stars  of  visual  Magnitude  *_  up  to  6.  Through 
a  telescope  of  2M  aperture  (a  likely  instrument  to  be  used  in  a 
satellite)  stars  of  visual  Magnitude  up  to  10  or  12  can  be  detected. 

The  wave  length  region  best  suited  for  ueasuresmrs  required  in 
satellite  navigation  is  discussed  in  Report  II.  It  is  approximately 
the  region  from  3000  A  to  5000  A.  On  the  other  hand  the  visual  asgni- 
tude  refers  to  stellar  brightness  in  the  spectral  region  of  greatest 
visual  sensitivity,  which  centers  at  5560  A  and  extends  fron  epprox- 
inately  5000  A  to  6200  A. 

Assuming  Planck's  law  to  hold  for  stars,  a  knowledge  of  the 
effective  surface  temperature  and  visual  Magnitude  can  be  used  to 
derive  total  irradiancc  ir.  any  other  spectral  region  (such  as  that  from 
3000  to  5000  A).  However,  for  this  prelialnary  discussion  this  is  not 
considered  necessary.  Instead  we  can  use  the  photographic  magnitudes 
as  a  criterion  of  brightness  in  the  above-mentioned  spect  ■»!  region, 
since  this  region  is  almost  coincident  with  the  regionof  maximum 
sensitivity  of  the  photographic  plate.  The  difference  between  photo¬ 
graphic  magnitude  ny  and  the  visual  magnitude  is  known  as  color 
index,  C.I. 


The  color  index  for  individual  bright  stars  is  given  by  Allen  (1957, 
Section  112),  and  for  various  spectral  and  luminosity  classes  it  is 
given  in  Table  1.  Right  ascension,  declination,  photographic  mag¬ 
nitudes  and  spectral  and  luminosity  classes  for  the  50  (photograph¬ 
ically)  brightest  stars  is  given  in  Table  3  (Arranged  in  decreasing 
order  of  brightness;  data  taken  from  Allen  1957,  Section  112) 
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TABLE  3 


FIFTY  BRIGHTEST  SUSS  (PBOIOGSAPuIC  MAGSHOTE) 


STAR 

a 

6 

■ 

PS 

Spec¬ 

tral 

Class 

Luminosity 

Class 

Sirius 

a  CM* 

6h40.8* 

-16*35* 

-1.58 

Al 

V 

Canopus 

a  Car 

6^21.7* 

-52*38* 

-0.67 

FO 

T 

Vega 

a  Lyr 

18h33.6- 

+38*41* 

-0.08 

AO 

V 

Rigel 

3  Ori 

5h  9.7* 

-  8*19* 

-0.03 

B8 

T 

Achernar 

a  Eri 

ln34.0* 

-57*45* 

0.27 

B5 

IV 

3  Cen 

13h56.8* 

-59*53* 

C.31 

B1 

II 

a  Cru 

12h21.0* 

-62*33* 

0.40 

B1 

IV 

a  Cen 

14h32.8® 

-60*25* 

0.43 

G2,  K3 

Procyon 

acxi 

7h34.l" 

+  5*29* 

0.66 

F5 

IV 

Spies 

a  Virgo 

I3h19.7,n 

-10*38 

0.68 

B1 

V 

Altair 

a  Aql 

19h45.9* 

+  8*36* 

0.87 

A7 

V 

Capella 

a  Aur 

5h  9.3” 

+45*54* 

0.87 

G1 

III 

3  Cru 

12h41.9n 

-59*  9* 

0.94 

BO 

III 

Segulus 

a  Leo 

10h  3.0ra 

+12*27’ 

1.10 

B8 

V 

Foraalhaut 

a  Ps  a 

22h52.l” 

-30“  9' 

1.13 

A3 

V 

£  C  Ha 

6h54.7C1 

-28*50’ 

1.20 

B2 

II 

Deneb 

a  Cyg 

2Ch18.6B 

+44*55’ 

1.20 

A2 

I 

Arcturus 

a  Bob 

14nll.Im 

+19*42' 

1.24 

KO 

III 

).  SCO 

i7h26.8m 

-37“  2’ 

1.27 

B2 

IV 

y  Ori 

5h19.8° 

+  6*16’ 

1.37 

B2 

III 

I 

i 
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TABLE  3 

(Continued) 

x» 

;  i 

i . 

STAR 

a 

-h _ ■ 

5 

a 

re 

Spec¬ 

tral 

Class 

Luminosity 

Class 

«  I 

V«»bVt  U  VC« 

/  zo.z 

+3a  6 ■ 

i.j» 

A1 

V 

g5  * 

e  Cri 

-  i°is* 

1  43 

B0 

I 

p  Tau 

5**20.0* 

+28=32’ 

1.46 

B8 

III 

»  i 

e  U  Ms 

12h49.6* 

+56*30* 

1.47 

Alp 

5  1 

r 

7  Vel 

8h  6.4* 

-47*  3’ 

1.51 

WG7 

i  I 

tjUHa 

13h43.6* 

+49*49’ 

1.56 

B3 

V 

£  1 

*. 

a  Gru 

22h  1.9* 

-47*27’ 

1.58 

B5 

V 

1  1 

■  p  Car 

9h12.1* 

-69*18’ 

1.59 

AO 

III 

'(  | 

5  Ori 

5h35.7* 

-  2“  O’ 

1.60 

09 

I 

i  1 

3  C  Ha 

6h18.3* 

-17*54’ 

1.63 

B1 

III 

| 

a  Pav 

20h17.7® 

-57“  3’ 

1.65 

B3 

IV 

€  Sgr 

18h17.5® 

-34*26’ 

1.75 

B9 

IV 

1 

a  Sgn 

18h49.l“ 

-26*25’ 

1.77 

B3 

V 

f 

6  Vel 

18h41.9® 

-54*21’ 

1.83 

A0 

V 

;  1 

7  Gem 

6h31.9® 

+16*29’ 

1.85 

A1 

IV 

i 

*  Ori 

5h43.0® 

-  9*42’ 

1.85 

B0 

I 

\ 

$  Pup 

8h  0.1® 

-39*43' 

1  86 

05 

i 

O  And 

0^  3.2® 

+28*32’ 

1.91. 

B8p 

III 

i 

5  Ori 

5h26.9® 

-  0*22’ 

1.91 

09 

II 

i 

P  Aur 

5h52.2® 

+44*56’ 

1.93 

A2 

IV 

i 

Algol  p  Per 

3h  1.6® 

+40*34’ 

1.95 

BB 

V 

1 

i 

i 

i 

Pollux  P  Gem 

7h39.2® 

+28*16' 

26 

2.00 

G«5 

III 

TABLE  3  (Continued) 


?. 

a 

5 

■ 

PS 

Spec¬ 

tral 

Class 

Lussix^sity 

Class 

K  Sco 

17"35.6= 

-38*59' 

2.00 

B2 

IV 

7  Cas 

0hjC.7* 

+50*11’ 

2.01 

BDe 

It 

e  Cen 

13h33.6* 

-52*57  * 

2.06 

B1 

IV 

7  Cen 

ttNe.o* 

-48*35’ 

2.08 

AO 

IV 

a  Cr  B 

15h30.4* 

+27“  3* 

2.10 

AO 

V 

3  Leo 

llh44.0° 

+15°  8’ 

2.12 

A3 

V 

i]  Cen 

14h29.2* 

-41*43’ 

2.13 

B3p, 

A2p 

5  3  Ma 

13h19.9® 

+55*27* 

2.14 

A2 

V 
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5.3  FREQUENCY  OF  STARS  OF  DIFFERENT  APPARENT  MAGNITUDES 

Figure  8  gives  the  number  of  stars,  Q  ,  brighter  than  a  given 

IS 

photographic  or  visual  Magnitude  (m^  or  u^) .  The  increase  of  0^ 
with  a  is  approximately  logarithmic,  but  tapers  off  for  higher  magni¬ 
tudes.  It  is  also  evident  froa  this  figure  that  the  number  of  stars 
brighter  than  a  given  photographic  magnitude  is  less  than  two-thirds 
the  number  brighter  than  the  same  visual  magnitude,  and  this  ratio 
decreases  as  vs  go  towards  fainter  stars,  becoming  less  than  one-naif 
at  magnitude  10,  (See  Table  4.) 

5.4  DISTRIBUTION  OF  STARS  OVER  THE  CELESTIAL  SPHERE 

Figure  9  st&vs  the  distribution  of  the  fifty  photographically 
brightest  stars  using  equatorial  coordinates.  It  is  evident  that  the 
distribution  is  far  from  uniform.  The  two  solid  lines  represent 
parallels  of  galactic  latitude  (+15°  and  -15°) .  Thirty-two  of  Che  fif¬ 
ty  stars  are  within  this  narrow  band. 


The  ron-uniformity  of  distribution  increases  rapidly  as  we  go 

towards  greater  and  greater  limiting  magnitudes,  as  is  shown  by  the 

plots  in  Figures  10  and  11  where  log  N  (number  of  stars  brighter 

in 

than  photographic  magnitude  ®  ,  per  sqBSre  degree)  is  plotted  against 

PR 

galactic  latitude  (for  m  *  5,  IQ,  15  and  20)  and  against  (for 
galactic  latitudes  10°,  ±30°,  ±60“  and  ±90°)  respectively.  It  is  evi¬ 
dent  that  for  any  limiting  photographic  magnitude  m  ,  there  are  fewer 

P8 

stars  at  high  galactic  latitudes  than  near  the  galactic  equator,  and 
that  this  ratio  decreases  with  increasing  limiting  magnitude. 


The  model  of  our  galaxy  as  a  lens-shaped  system  with  the  sun 
located  in  an  eccentric  position  near  the  galactic  plane  explains  the 
observed  concentration  of  stars  in  low  galactic  latitudes. 
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Figure  8 
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magnitude  m— *• 


The  number  of  stars  brighter  then  apparent 
magnitude  m. 
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TABLE  4 

NUMBER  OF  STARS  BRIGHTER  THAN  APPARENT  MAGNITUDE  tn 


(VISUAL  AND  PHOTOGRAPHIC) 

tn 

0 

n 

m 

m 

(photographic) 

(visual) 

4 

3.236  x 

102 

5.370  x  10: 

6 

3.020  x 

103 

4.898  x  10 

8 

2.344  x 

104 

4.169  x  10J 

10 

1.738  x 

10s 

3.467  x  10; 

12 

1.202  x 

106 

2.455  x  101 

14 

7.586  x 

106 

1.514  x  10 

16 

3.981  x 

107 

8.318  x  10 

18 

1.820  x 

108 

3.802  x  101 

20 

6.166  x 

108 

1.175  :<  101 

30 
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Figure  11,  Number  of  stars  brighter  chan  magnitude  m  for  various 


values  of  m. 


I 


I 


\ 


SECTION  6 

FREQUENCY  OF  STAR  OCCULTATIONS 

6.1  THE  FRACTIONAL  AREA  OF  A  CELESTIAL  SPHERE  OCCULTED  INSTANTANEOUSLY 

The  semi-sngle  of  the  earth's  "shadow  cone",  according  to 
Equation  (3.1),  is 


This  corresponds  to  a  solid  angle 
w  *  2«(1  -  cos  $) 


(6.1) 


The  fractional  area  of  a  celestial  sphere  occulted 
-  ^  *  cos  (3) 

"  \ 


r  a2  1 

l  - 

1  -  — 

2 

1 

r 

/ 

h 


1  - 


1  - 


1  - 


a2(l  -  e)2J 


a 


1  - 


a2(l  +  e)2 


at  the  perigee 


at  the  apogee 


\  L 

6.2  THE  FRACTIONAL  AREA  OF  A  CELESTIAL  SPHERE  OCCULTED  PER  PERIOD 


Consider  the  celestial  sphere  drawn  in  Figure  12  where  the  cen¬ 
tral  plane  is  that  of  the  satellite  orbit.  The  band  of  the  celestial 
sphere  enclosed  between  the  dotted  parallels  of  latitude  denotes  the 
shadow  band  for  a  circular  orbit.  However  for  an  elliptical  orbit  the 
band  will  not  have  a  uniform  thickness.  The  maximum  thickness  will 
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Figure  12.  Area  of  celestial  sphere  occulted.  The  thickness  of  the 
"shadow  band"  varies  during  an  orbital  revolution  for  an 
elliptical  orbit. 
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occur  ac  the  apogee  psinr  fccrresposdisg  to  the  satellite  location  at 
the  perigee)  and  the  wniwe  at  the  perigee  point.  The  shaded  area 
represents  such  a  shadow  band  c£  varying  thickness.  It  is  drawn  to 
correspond  to  an  ellipse  whose  perigee  distance  equals  the  radius  of 
the  circular  orbit  whose  ’’shadow  cone”  is  represented  by  the  dotted 
parallels  of  latitude.  The  fractional  area  occulted  per  period  is 
given  by: 


4s*2  ^o 


f+e 


r  *2  sin  S  dS  dv 


u[-~  f 

2  " 


i-  r' 


2  sin  pdv 


^  p~  ao(l  -  e  cos  v) 

*J_  °  .(1  - 


i  ao  r 

2"  a(l  -  e2)  L 


e  sin  v 


a(l  -  e4) 


i  a 
1  _ o 

1  +  e  r 
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where  r  *  s(l  -  e)  i*  the  perigee  distance.  Toe  average  fractional 
P 

area  occulted  per  boor  is. 


(6.4) 


where  P  is  given  by  Equation  (2.3). 
o 

For  a  circular  orbit  of  radius  r^,  the  average  fractional  area 
occulted  per  hour  is: 


1  +  e 

which  is  - ^rjz  tiroes  as  large  as  tor  an  ellipse  whose  perigee 

(1  -  *)' 

distance  equals  the  radius  of  the  given  circle.  For  c  =  ^  this 
factor  equals  4.2 

It  is  to  be  noted  that  there  can  be  a  fairly  substantial  variation 
in  the  actual  rate  of  occultation  of  the  fractional  area  of  the  celes¬ 
tial  sphere  on  account  of  the  fact  that  the  speed  of  the  satellite 
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as  mil  as  the  instantaneous  fractional  area  occulted,  have  their 
aodam  values  ac  res  perigee  sac  their  slntsan  values  at  the  apogee. 

6.3  OCO-LTATIOT  FREOUESCY 

We  saw  in  Section  4,  that  the  ainiaa  njafcer  of  occultation 
observations  of  distinct  stars  required  rc  calculate  all  six  orbital 
elements  is  six.  Since  experiaental  errors  are  inherent  to  all  sscas- 
ureaents,  it  would  be  desirable,  in  order  to  reduce  the  probable 
errors  c£  Measurements,  to  make  approxiaately  twelve  observations  and 
obtain  least  square  solutions  for  the  unknown  elements.  This  seens 
particularly  desirable  in  view  of  the  variable  and  somewhat  unknown 
radiance  (due  to  aurorae,  air  glow,  etc.)  of  the  earth’s  atmosphere 
through  which  the  star  light  would  have  to  pass  on  its  way  to  the 
satellite  measuring  device.  This  particular  problem  is  discussed  in 
another  report  of  this  scries  (Jones  and  Kaqvi  1962,  hereafter  refer¬ 
red  to  as  Report  111). 

We  expect  that  in  most  cases,  the  approximate  orbital  elements 
would  be  known  even  before  the  launch,  and  this,  in  conjunction  with 
the  condition  of  occultation  derived  in  Section  3,  would  permit  us  to 
calculate  exactly  which  (and  how  many)  stars  will  be  available  for 
occultation  observations.  However,  some  general  results  giving  the 
frequency  of  occultation  can  be  obtained  if  we  assume  a  uniform 
distribution  of  stars  over  the  sky.  The  fact  that  the  stars  are 
concentrated  in  low  galactic  latitudes  (See  Section  5)  will  result, 
depending  upon  the  inclination  of  the  orbital  plane,  in  a  higher  or 
lower  occultation  frequency  than  that  derived  on  the  basis  of  an 
assumed  uniform  distribution  of  the  same  number  of  stars.  If  the 
orbital  plane  is  tnclined  at  a  small  angle  to  the  galactic  plane, 
more  occultations  would  occur  than  predicted,  and  vice  jersa.  Since 
it  would  be  desirable  to  have  twelve  occultations  even  if  the  orbital 
plane  were  at  right  angles  to  the  galactic  plane,  it  is  estimated 
that  a  total  of  approximately  25  star  occultations  should  be  allowed 
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for  if  we  assuae  a  unifora  distribution  of  stars.  Very  approximately 
it  would  scan  that  for  certain  favorable  orbits  (low  inclinations  to  the 
galactic  plane)  the  fulfillment  of  the  conditions  would  give  us  fifty 
or  aore  stars  whereas  for  the  least  favorable  orbits  there  will  be 
about  twelve  stars  available. 

Let  us  assuae  that  H  stars  are  uniformly  distributed  over  the 
sky.  Froa  Equation  (6.3),  the  niaber  of  stars  occulted  per  period  of 
revolution  is 


f 

P 


N 


1  +  e  r 


(6.6) 


The  average  number  of  stars  occulted  per  hour,  according  to 
Equation  (6.4)  is 


N_ 

P 

o 


U  -  e) 

1  +  e 


3/2 


(6.7) 


In  order  that  q  seer  occultations  can  be  observed  per  satellite 
period,  the  total  number  of  stars  available  over  the  entire  sky 
should  be 


r 

N  *  q(l  +  e)  — ^ 
«  ao 


(6.8j 


Similarly  if  we  »-'58h  to  observe  an  average  of  q  star  occultations 
per  hour,  the  total  number  of  stars  over  the  entire  sky  should  be 
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01004-1250? 


Figure  13.  Number  of  stars  required  and  the  corresponding  apparent 
photographic  magnJtudes  for  25  occultations  per  period 

'it 

(Njj)  and  per  hour  (N  for  a  circular  orbit  (e  «=  0). 
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APPARENT  PHOTOGRAPHIC  MAGNITUDE 


1  +  e 


(6.9) 


* 

K 


q 


q? 


(1  -  e) 


3/2 


In  Figure  13.  we  have  plotted  and  N*  respectively  against 

h  (*  r  -  a  )  for  circular  orbits  (e  «  o)  and  for  q  *  25.  For  orbits 
p  p  o 

of  any  other  eccentricity  these  values  should  be  Multiplied  with  (1  +  e) 
1  -t-  e 

and  r - ryr  respectively.  The  Magnitude  scale  corresponding  to 

(1  -  e) 

the  nuMber  of  stars  is  given  on  the  right  side  of  the  graph. 


Although  the  present  report  deals  prinarily  with  the  case  neglect¬ 
ing  the  atmosphere,  there  are  two  effects  due  to  the  atmosphere, 
discussed  in  Reports  II  and  III,  which  affect  the  occultation  frequency 
and  the  nunber  of  stars  required  (or  the  limiting  stellar  magnitude), 
which  should  be  mentioned  here.  First,  due  to  Rayleigh  scattering, 
the  sunlit  half  of  the  earth's  atmosphere  has  such  a  large  intensity 
that  no  useful  observations  of  stellar  occultation  against  this  bright 
foreground  are  possible.  As  a  result  the  frequency  of  occultations 
is  effectively  reduced  by  ~  .  The  minimum  requirement  on  the  number 
of  stars  is  increased  by  2,  and  the  limiting  magnitude  of  the  detecting 
instrument  by  approximately  0.65  magnitude.  Second,  due  to  partial 
extinction  of  star  light  (again  due  to  Rayleigh  scattering)  the  light 
intensity  is  decreased  by  approximately  a  factor  of  two,  at  the 
instant  of  observation  (See  Report  II).  This  does  not  affect  the 
frequency  of  occultations,  or  the  requirement  on  the  number  of  stars, 
but  does  increase  the  limiting  apparent  magnitude  for  which  the  instru¬ 
ment  should  be  designed  by  approximately  0.75.  When  the  two  effects 
are  combined,  the  requirement  on  limiting  magnitude  is  increased  by  1.4. 
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TABLE  5 


NUMBER  OF  STARS  REQUIRED  AND  THE  CORRESPONDING  INSTRUMENTAL  LIMITING 
MAGNITUDE  PCS  25  QCCULTATIONS  PER  PERIOD  <H25)  AND  PER  HOUR  (N*_) 

FOR  A  FEW  SELECTED  ORBITS 


Perigee  Height  h  *  200  ion 


Eccen¬ 

tricity 

semi- 
major 
axis 
a  (ks) 

Period 

(hours) 

Number 

of 

Stars 

N25 

Limiting 

Magnitude 

(photo¬ 

graphic) 

* 

N25 

Limiting 

Magnitude 

(photo¬ 

graphic) 

0 

6578 

1.475 

26 

1.58 

38 

1.91 

.01 

6644 

1.497 

26 

1.58 

39 

1.93 

.5 

13156 

4.171 

39 

1.93 

161 

3.38 

Parigee 

Height  h 

*  2000  km 

0 

8378 

2.106 

33 

1.83 

70 

2.55 

.01 

8463 

2,152 

33 

1.83 

72 

2.56 

.5 

16756 

5.995 

49 

2.14 

295 

3.92 

NOTE:  (I)  Circular  and  elliptical  orbits  are  considered  in  this 
table  but  only  the  circular  orbits  arc  considered  in 
Figure  13. 

(2)  The  limiting  magnitude  may  be  raised  by  as  much  as  1.4 
magnitudes  on  account  o£  extinction  and  radiance  of  the 
earth's  atmosphere  (See  end  of  this  section). 
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